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Abstract 

We investigate the allowed configurations in the stationary state of the cellular au- 
tomaton model for single-lane traffic. It is found that certain states in the configura- 
tion space can not be reached if one uses parallel dynamics. These so-called Garden of 
Eden (GoE) states do not exist for random-sequential dynamics and are responsible 
for the strong short-ranged correlations found in parallel dynamics. By eliminating 
the GoE states we obtain a simple and effective approximative description of the 
model. For v max = 1 the exact solution is recovered. For v max = 2 this elimination 
leads to much higher values of the flux compared to the mean-field result which are 
in good agreement with Monte Carlo simulations. 
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1 Introduction 



The description of traffic flow using cellular automata (CA) is quite successful 0, 
despite the simplicity of the model 0. CA are, by design, ideal for large-scale 
computer simulations. This fact has already been used for the simulation of urban 
traffic in various cities, see e.g. |3], f|, [5]. On the other hand, analytical descriptions are 
difficult. In || |7|, || we have developed several methods which yield an approximate 
description of the stationary state. In certain limits (e.g. v max = 1 or p — > 0), these 
methods even become exact (see || for a review). These approaches are based on 
a microscopic description which takes into account certain correlations. Part of the 
difficulties come from the fact that one uses parallel dynamics. This introduces a 
non-local aspect into the problem since the whole lattice is updated at once. On the 
other hand, random-sequential dynamics is much simpler to treat analytically. For 
Vmax = 1 for instance, simple mean-field theory gives already the correct steady state, 
i.e. there are no correlations. Here we propose a rather simple analytical approach 
which exhibits the main difference between parallel and random- sequential dynamics 
very clearly. For parallel dynamics not all states of the configuration space can be 
reached by the dynamics, some are 'dynamically forbidden'. This is not the case for 
random- sequential dynamics. 

For completeness we recall the definition of the CA model for single-lane traffic ||. 
The street is divided into L cells which can be occupied by at most one car or be 
empty. The state of each car is described by an internal parameter ('velocity') which 
can take on only integer values v = 0,1,2, ... , v max . The state of the system at time 
t + 1 can be obtained from the state at the previous time t by applying the following 
four simple update steps to all cars at the same time (parallel dynamics): 

Rl Acceleration: Vj(t) — > Vj(t + |) = mm{vj(t) + 1, v max } 

R2 Braking: if Vj(t+±) > dj(t) then Vj(t + |) -> Vj(t + §) = dj(t) 

R3 Randomization: Vj(t + §)—>■ Vj(t + 1) = max{0, Vj(t + |) — 1} 

with probability p 

R4 Driving: car j moves Vj(t + 1) cells. 

Here dj(t) denotes the number of empty cells in front of car j, i.e. the gap or headway. 
In the following it will be important that Vj(t) is just the number of cells that car j 
moved in the timestep t — 1 — > t. 

In we investigated a microscopic mean-field theory (MFT) for this kind of cellular 
automaton. The most important result was that the flows obtained are much too 
small compared to computer simulations. The reason is that MFT cannot account 
for the 'particle-hole attraction' found in the stationary state, i.e. the probability to 
find an empty cell in front of a (moving) car is enhanced compared to the mean-field 
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result. This effect is taken into account by the n— cluster approximation introduced 
in H 0| . Here the lattice is divided into clusters of length n which overlap n — 1 cellsQ. 
It turned out that the 2-cluster approximation is exact for v max = 1. For v max = 2 
already small cluster sizes (e.g. n = 4) gave an excellent agreement with numerical 
results for the flow. 

In H an alternative analytic approach was introduced, the so-called car-oriented 
mean-field theory (COMF). Here the dynamical variables are not the occupancies of 
the cells but the gaps dj(t) between consecutive cars. In the COMF these gaps are 
treated as being independent. Again we found that COMF is exact for v max = 1. 
For v m ax = 2 and the 'traditional' value p = 0.5 the flow obtained from COMF is 
comparable to the 3-cluster result. However, the COMF seems to become exact in 
the limit p — > 0, in contrast to the 2-cluster approach ||. 

In the present paper we like to present a rather simple extension of MFT. The key 
idea is a reduction of the configuration space by removing all states which cannot 
by reached dynamically. In the context of cellular automata these states are called 
Garden of Eden ( GoE) states or paradisical states. A simple example (see Fig. |l|) for 
Vmax = 2 is the configuration (•, 1,2) of two consecutive cells, where '•' denotes an 
empty cell and the numbers correspond to the velocities of the cars. Cars move from 
left to right. Obviously the velocity is just the number of cells the car moved in the 
previous timestep. Therefore the configuration (•,1,2) could have evolved only from 
a state which has two cars in the leftmost cell. Since double-occupations are not 
allowed in the present model, states containing (•,1,2) are dynamically forbidden, 
i.e. are GoE states. 
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Figure 1: A Garden of Eden state for the model with v max > 2. 

Here we reinvestigate MFT for v ma x = 1 and v ma x = 2, but eliminate the GoE states. 
This paradisical mean-field (pMF) theory will lead to a considerable improvement of 
the results. 

2 Mean-field theory 

Here we briefly review the MFT results for v ma x = 1 and v ma x = 2. A more complete 
account along with a detailed derivation can be found in 0. The densities of cars 

^MFT correponds to n = 1. 
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with velocity v (v — 0, 1, . . . , v max ) in the stationary state is denoted by c v . Therefore 
the full density of cars is given by c = YZ=o x °v- F° r convenience we also introduce 
the 'hole' density d = 1 — c and the abbreviation q = 1 — p. The flux (or current) is 
given by /(c) = £"=f vc v . 

For v m ax = 1 the mean-field equations then read: 

c = (c + pd)c, (1) 
ci = gcrf (2) 

so that the flux is simply given by 

/SUc) = Ci = gc(l-c). (3) 



For v max = 2 the rate equations for the densities are given by 



The solution 



c = (c + pd)c + (1 + pd)c(c 1 + c 2 ), (4) 
Ci = d [qc + (qc + pd) (d + c 2 )] , (5) 
c 2 = gd 2 (ci + c 2 ). (6) 



_ (l+^)c 2 



Cl 



g(l — qd 2 )dc 
1 — p<i 2 



yields for the flux 



c 2 = ^ (9) 



ftikc) = c 1 + 2c 2 = ?^J^. (10) 



Comparison with the results from Monte Carlo simulations shows that the MFT 
results underestimate the flow considerably 0. Therefore strong short-range corre- 
lations exist which increase the flow compared to the prediction of MFT. Using the 
exact solution for the case v max — 1 @ one can demonstrate explicitely the existance 
of a strong particle-hole attraction, i.e. the probability to find an empty cell in front 
of an occupied cell is enhanced compared to the MFT result. 

At first sight, this result is surprising. For random-sequential dynamics^ it is known 
that MFT is exact for v max = 1. Therefore the origin of the correlations is the 
parallel update procedure. In the following we will see that the existence of the 
GoE states is responsible for the differences between parallel and random-sequential 
dynamics. 

Hn random-sequential dynamics in each timestep a cell which is updated is picked at random. 



4 



GoE states for v max = 1 



The question, whether a state is a GoE state or not, can be decided locally by 
investigating just nearest-neighbour configurations. It turns out that GoE states are 
all states containing the local configurations (0, 1) or (1, 1), i.e. a moving vehicle is 
directly followed by another car. This is not possible as can be seen by looking at the 
possible configurations at the previous timestep. The momentary velocity gives the 
number of cells that the car moved in the previous timestep. In both configurations 
the first car moved one cell. Therefore it is immediately clear that (0, 1) is a GoE 
state since otherwise there would have been a doubly-occupied cell before the last 
timestep. The configuration (1, 1) is also not possible since both cars must have 
occupied neighbouring cells before the last timestep too. Therefore, according to 
rule R2, the second car could not move. 

We now modify the MFT equations (|I]) and (0) to take into account the existence of 
GoE states. Following the procedure described in J7| it turns out that only (P has 
to be modified. Due to this modification the normalization cq + C\ = c is no longer 
satisfied automatically. Therefore a normalization constant M has to be introduced. 
The final equations are than given by 

c = M{c +pd)c, (11) 
Cl = Nqcd, (12) 

with the normalization 

N = — (13) 
c + a 

Since Co + Ci = c we have only one independent variable for fixed density c, e.g. c\. 
Solving (|ITD, ( p~2"D for c\ we obtain 



d = i (l - ^/l-4 g (l-c)c) . (14) 

The flow is given by f(c) = c\ and we recover the exact solution for the case v max = 1 
§• 

This result confirms the expectations mentioned above. One can see clearly that the 
difference between random-sequential and parallel dynamics is the existence of GoE 
states in the latter. After eliminating these GoE states, no correlations are left in 
the reduced configuration space. 



4 GoE states for v max = 2 

In this case more GoE states exist. In order to identify the GoE states it is helpful 
to note that the rules R1-R4 imply dj(t) — dj(t — 1) + Vj+i(t) — Vj(t) and therefore 

dj{t) > v ]+1 {t)- Vj {t), (15) 
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Vj (t) < djit-l) (16) 
The second inequality ([H]) is a consequence of R2. 

In the following we list the elementary GoE states, i.e. the local configurations which 
are dynamically forbidden (cars move from left to right): 

(0,1), (0,2), (1,2), (0,.,2), (17) 
(1,1), (2,1), (2,2), (1,.,2), (2,.,2), (18) 

(0,.,.,2). (19) 

The numbers give the velocity of a vehicle in an occupied cell and • denotes an 
empty cell. Using Monte Carlo simulations we have checked that there are no further 
elementary GoE states for clusters up to 10 cells. 

The elementary GoE states in (|i~7|) violate the inequality fll~5|), and the configurations 
in (|T8| ) violate (|l6l). The state in (|19"D is a second order GoE state. Going one step 
back in time leads to a first order GoE state since (0, •, •, 2) must have evolved from 
(0, v) (with v = 1 or v = 2). 

Again we can derive the pMF equations by modifying the method for the derivation 
of the MF mean- field theory |7j. Taking into account only the first order GoE states 
(|TT|) and flTB| ) one obtains the following pMF equations: 

c = Af[c c + pd(c + cic)] , (20) 
ci = N pd 2 (ci + c 2 ) + qd(c + c x c) , (21) 
c 2 = Mqd 2 ( Cl + c 2 ). (22) 

The normalization M ensures Cq + c\ + c 2 = c and is given explicitly by 

j\f = — (23) 

c + dc\ + d 2 c 2 c + d{l — c 2 ) ' 

Using (22) we can express c 2 through Co and c only: 

C2 = Yd ( C ° + d ~ \/( c o + rf ) 2 - 4 ^ 3 ( c - c o)) • ( 24 ) 

Inserting this result into ( 2D ) we obtain a cubic equation which determines Co in 
terms of the parameters c and p: 

acl + pc 2 + 7C + 5 = (25) 

where the coefficients are given by 

a = —qd? — pd 2 + qd — q, (26) 

= (3p 2 - p - l)d A + (-p 2 -3p+ l)d 3 + (-p 2 +p + 2)d 2 -3qd + q, (27) 

7 = pc d[(-p 2 -3p + 2)d 3 + (p 2 +p+l)d 2 + (2p-3)d + q], (28) 

S = -(pedf. (29) 
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In principle the zeroes of (|25| ) can be determined explicitely. The physical solution is 
the one in the interval [0, 1]. The flow can be calculated as f{c,p) = c\ + 2ci where 
ci and C2 are determined by (^4j) and c = Co + c\ + C2. 




Figure 2: Fundamental diagram for i> maa; = 2 and p = 0.5 (left) and p = 0.1 (right). 
Comparison of paradisical MF (full line) with results from computer simulations (•) 
and the naive MF approximation (dotted line). 

We have also calculated the pMF equations which in addition take into account the 
second order GoE state fllTf). One has to modify the equations for cq and c\ only (and 
the normalization) by replacing the last term in ( ^0|) by pd{(co + ci)c + cod(l — C2)} 
and the last term in (|21|) by gd{(c + Ci)c + c <i(l — c 2 )}. However, these modifications 
lead only to minor changes in the results for the fundamental diagram. 



5 Discussion 

We have presented an analysis of the allowed configurations in the CA model for 
traffic flow. Due to the use of parallel dynamics not all configurations can be reached 
through the dynamics. Eliminating these Garden of Eden states allowed us to im- 
prove the results of the naive mean-field theory considerably. 

GoE states can be characterized locally. We identified all elementary GoE configu- 
rations for v max = 1 and v max = 2. It turns out that for v max = 1 it is sufficient to 
investigate only configurations of all clusters of two cells. For v max = 2 the largest 
elementary GoE configuration consists of four neighbouring cells. 

For v max = 1 the paradisical mean-field theory is able to reproduce the exact solution. 
This implies that in the subspace without GoE states all configurations are equally 
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probable. This has to be compared with random-sequential dynamics. Here all 
configurations are equally probable and naive mean-field theory is exact. This means 
that the strong short-ranged correlations found for parallel update are solely due to 
the use of parallel dynamics. 

In fact one may speculate that this is rather general. The difference between random- 
sequential and parallel update comes mainly from the existence of GoE states in the 
latter. This implies that a method that 'works' for random-sequential dynamics (e.g. 
an exact solution or good approximation) should also work for parallel dynamics, but 
now in the subspace without GoE states. 

For v max = 2 the paradisical mean-field theory yields a considerable improvement 
of the mean-field results, but it does not become exact. One observes a qualitative 
difference to the case v max = 1, since now there are correlations present which can 
not be explained by the existence of GoE states. 

The existence of GoE states gives a simple criterion for the quality of an approxi- 
mation: A good approximation should be able to account for all GoE states. This 
can be illustrated for the case v max = 1. The methods used previously for the ex- 
act solution are the 2-cluster approach || |7| and Car-Oriented Mean-Field theory 
(COMF) [§. Both methods are able to identify both GoE states (0,1) and (1,1). For 
Vmax = 2 one needs at least the 4-cluster approximation to account for all GoE states 
( p!7| - [T9|) . Indeed, the results of || [jj show that the 4-cluster results are in excellent 



agreement with simulation results. COMF is able to identify all GoE states since all 
elementary GoE configurations consist of only two neighbouring cars, i.e. there are 
no elementary GoE states with three vehicles. 

Finally we want to point out that the model is ergodic in the sense that for configu- 
rations r, t' which appear in the stationary state with probabilities P(r), P(r') > 
there is a non- vanishing transition probability P(t — > r 1 ). The existence of GoE 
states poses no problems in computer simulations. If the initial state is a GoE state, 
it will become a non-GoE state after the first timestep. 

The method presented here is also applicable to other models. An interesting case is 
the asymmetric exclusion process (see |H| and references therein) , which is identically 



to the model investigated here with v max = 1, but with open boundary conditions 
where an injection and/or removal of particles is possible. In [[□]] the existence of 
GoE states has been used to obtain an approximative description of the deterministic 
limit of the ASEP which is in excellent agreement with numerical results. 
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